Meron- Cluster Simulation of the Vacuum in 



W. Bietenholz , A. Pochinsky a ' & and U.-J. Wiese a 

a Center for Theoretical Physics, 
Laboratory for Nuclear Science, and Department of Physics 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139, U.S.A. 

h Institute for Theoretical and Experimental Physics 
117259, Moscow B. Cheremushckinskaya 25, Russia 



The 2-d 0(3)-model with a ^-vacuum term is formulated in terms of 
Wolff clusters. Each cluster carries a half-integer topological charge. The 
clusters with charge ±1/2 are identified as merons. At 6 = tt the merons 
are bound in pairs inducing a second order phase transition at which 
the mass-gap vanishes. The construction of an improved estimator for 
the topological charge distribution makes numerical simulations of the 
phase transition feasible. The measured critical exponents agree with 
those of the k = 1 Wess-Zumino-Novikov-Witten model. Our results are 
consistent with Haldane's conjecture for 1-d antiferromagnetic quantum 
spin chains. 
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Some time ago Haldane conjectured [1] that integer and half-odd-integer 1-d an- 
tiferromagnetic quantum spin chains behave qualitatively differently. While integer 
spin chains have a mass-gap half-odd-integer chains should be gapless. This has 
been confirmed numerically for finite chains of spin f and 2 [2, 8] and analytically 
for half-odd-integer spins and for spin f [4]. The long-range physics of f-d quantum 
spin chains is described by an effective 2-d classical 0(3)-model. Haldane argued 
that the effective action for a chain of spins S contains a topological term i6Q. Here 
Q is the topological charge and 6 = 2ttS is the vacuum angle. Since the physics is 
periodic in 6, i.e. 6 G] — 7r,7r], integer spins have 6 = and half-odd-integer spins 
have 6 = tt . The standard 0(3)-model has a mass-gap in agreement with Haldane's 
conjecture. On the other hand, the conjecture implies that the mass-gap of the 
0(3)-model disappears at 6 = tt. This corresponds to a phase transition in the vac- 
uum angle. Due to the complex action it is notoriously difficult to simulate #-vacua 
numerically. A previous numerical study that was limited to \6\ < 0.8 tt found no 
phase transition in that region [5]. In fact, Haldane's conjecture has not yet been 
verified in the context of the 0(3)-model. In this paper we use the Wolff cluster 
algorithm [6] combined with a reweighting technique [7] to attack this problem. The 
construction of an improved estimator for the topological charge distribution enables 
us to simulate 6- vacua reliably for any value of 6. 

Affleck and Haldane have suggested a dynamical mechanism that explains why 
the mass-gap disappears at 6 = tt [8]. In this picture pseudo-particles with topologi- 
cal charge ±1/2 — so-called merons — are the relevant degrees of freedom. At 6 = 
the merons form an ideal gas. They disorder the system and thereby give non-zero 
mass to the physical particles. At 6 = 7r, on the other hand, the merons are bound 
in pairs and thus do not generate mass. Affleck confirmed this picture in a model 
where the 0(3) symmetry is explicitly broken to 0(2). Then the merons behave 
like vortices and the phase transition in 6 is analogous to the Kosterlitz-Thouless 
transition of the 0(2)-model. When the explicit 0(3) breaking is switched off it 
is unclear if this dynamical picture still holds. In fact there exists no definition of 
merons beyond the semiclassical approximation. 

Here we formulate the 0(3)-model in terms of Wolff clusters. Using an appro- 
priate action on a triangular lattice each cluster has a uniquely defined half-integer 
topological charge. It turns out that most clusters are neutral, some have charges 
±1/2 and very few carry larger charges. It is natural to identify merons as Wolff 
clusters with charge ±1/2. Indeed at 6 = the Wolff clusters are completely inde- 
pendent, as they should in order to resemble merons. At 6 = 7r, on the other hand, 
no clusters with half-odd-integer charges persist. Instead they form bound pairs of 
integer charge. Hence the cluster formulation of the model provides a definition of 
merons beyond the semiclassical approximation, and it confirms Affleck's dynamical 
mechanism in the 0(3)-model even without introducing symmetry breaking terms. 

Actually Affleck's picture of the phase transition is more quantitative [9]. He 
argues that the model at 6 = tt is a conformal field theory in the same universality 
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class as the k = 1 Wess-Zumino-Novikov-Witten model [10]. The critical exponents 
of this model are known analytically. The results of our numerical simulations are 
consistent with them, confirming Haldane's conjecture as well as Affleck's picture in 
the context of the 0(3)-model. 

For technical reasons we work on a 2-d triangular lattice with V = 3L 2 points. 
The lattice covers the volume V of a regular hexagon with periodic boundary con- 
ditions. To each lattice site x we attach a classical 3-component spin e x £ S 2 (of 
unit length). The action 

= £.(«„?.), <«,«)-{ jf"*^ ^1 >_1/2 ID 

{xy) ^ 

is a sum over nearest neighbor pairs (xy). It constrains neighboring spins to rela- 
tive angles less than 2tt/3 which is necessary for our purposes as will be explained 
below. The spin field defines a mapping from the lattice to the sphere S 2 . Using an 
appropriate interpolation, these mappings fall into topological classes characterized 
by an integer valued topological charge Q £ Z. The charge Q counts how many 
times S 2 is covered by the mapping. The total charge 

Q[e] = q{ex,e y ,e z ) , q(e x ,e y ,e z ) = A(e x , e y , e z ) /Air , (2) 

(xyz) 

is a sum of contributions from oriented elementary lattice triangles (xyz). The 
spins e X} e y and e z at the three corners x, y and z of a lattice triangle define a 
minimal spherical triangle on the sphere S 2 . Each lattice triangle contributes the 
area A(e X} e y} e z ) of the corresponding spherical triangle divided by the area of S 2 
with a sign depending on its orientation relative to (xyz). Another quantity of 
interest is the total magnetization M[e] = J2x^x- The partition function of a 6- 
vacuum is given by 

Z(6) = j Ve exp(-5'[e] + i6Q[e}) = Il f s2 d ^ exp(-5'[e] + i6Q[e}). (3) 

We also consider the topological susceptibility Xt(^) = ((Q 2 )e — (QYq)IV an d the 
magnetic susceptibility Xm(@) = ((M 2 )e — (M) 2 )/V. Using finite size scaling one 
can extract critical exponents from Xt(^) an d Xm(@)- 

The Wolff cluster algorithm is a very powerful tool that has been used to beat 
critical slowing down in numerical simulations of 0(N)-models [6]. Each step of the 
algorithm begins with the selection of a random direction characterized by a unit 
vector r. Then the spin components in this direction are updated by forming spin 
clusters which can only be flipped collectively. When flipped, a spin e x changes to 
e*E ' = e x — 2(r • e x )r. To define the clusters, nearest neighbor spins e x and e y are 
connected by a bond with probability 

p = 1 - min[l,exp(s(e :E , e y ) - s(e x ', e y ))]. (4) 
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With our definition of the action in eq.(l) a bond is always put if a spin flip would 
increase the relative angle of two neighboring spins beyond 2tt/3. Spins which are 
connected by bonds belong to one cluster. In the multi-cluster algorithm all spins of a 
cluster are flipped simultaneously with probability 1/2. The clusters are completely 
independent. This allows the construction of improved estimators for quantities like 

Xm(0) [6]. 

Here we use the cluster algorithm to update the ^-vacuum of the 0(3)-model. 
Actually we update at 6 = and include the extra Boltzmann factor exp(iOQ) in the 
measured observables. This requires a very precise determination of the topological 
charge distribution 

P(Q) = J V ^ ^Q,Q[e] exp(-5'[e]). (5) 

The partition function is then given by Z(Q) = J2q p(Q) exp(iOQ). Observables are 
computed separately in each topological sector 

(O)q = P(Q)- 1 J Ve S Q , Q[gi 0[e] exp(-S[e}), (6) 
and their ^-vacuum expectation value is given by 

(0) e = Z{e)- 1 Y.P{Q){0) Q exp(zeQ). (7) 

Q 

The distribution p(Q) varies over many orders of magnitude. In particular, large 
charges are suppressed exponentially. Still, at large values of 6 the contributions 
from the large charges are important and need to be determined precisely. In this 
respect a reweighting technique [7] using a trial distribution pt(Q) has been useful. 
The trial distribution should be as close as possible to the real distribution p(Q). In 
general this is difficult to achieve because one has not yet determined p(Q). In our 
case the improved estimator described below provides a good trial distribution pt(Q)- 
Then one works with an effective action S'efrfe] = S[e] + lnp t (Q[e]) and determines 
p(Q) = Pt(Q) I Sq,q[?] exp( — S'efrfe ]). This method can be combined with the 
multi-cluster algorithm. The clusters are grown as usual. Flipping a specific cluster 
changes the topological charge from Q to Q' . Hence we define the cluster charge as 
(Q — Q')/2. Therefore each cluster carries a half-integer topological charge. Due to 
the triangular lattice and our choice of action — which constrains the relative angles 
of neighboring spins — the cluster charges are uniquely defined, i.e. the charge of a 
specific cluster is independent of the orientations of all other clusters. Moreover the 
cluster charges can be computed locally. These are highly nontrivial consequences 
of our construction. The trial distribution pt(Q) enters only when the clusters are 
flipped. The probability to flip a specific cluster is given by the Metropolis algorithm 

J mm[l,p t (Q)/ Pt (Q')] if p t (Q>) + p t {Q), 
1 1/2 otherwise. 

Since the clusters are independent we can construct an improved estimator for the 
topological charge distribution. This is the essential ingredient of our method to 
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simulate a ^-vacuum. A configuration with n clusters represents an ensemble of 2 n 
configurations related to the original one by cluster flips. The weight of a configura- 
tion depends on its topological charge and is given by Pt{Q)~ X ■ Knowing all cluster 
charges, it is a simple combinatoric exercise to construct the topological charge dis- 
tribution for the ensemble of configurations analytically. This increases the statistics 
by a factor 2 n . Typically there are 0(100) clusters in our configurations which re- 
sults in a tremendous improvement. In particular, the partition function Z(tt) gets 
a positive contribution from each configuration. Hence there is no sign-problem un- 
like in other numerical methods. Similarly, we construct an improved estimator for 
the magnetic susceptibility. This is done separately for each topological sector. The 
orientation of the charged clusters determines the topological charge Q. When the 
neutral clusters are flipped we stay in the given topological sector. The improved 
estimator for (M 2 )q includes all configurations that result from flipping the neutral 
clusters only, and it is otherwise identical with the standard improved estimator [6]. 
It turns out that most clusters are neutral and the improvement factor is still large. 

We emphasize that the cluster formulation of the 0(3)-model provides more 
than just an efficient numerical algorithm. The clusters represent physical objects 
that play a crucial role in the dynamics. Indeed we identify the clusters of charge 
±1/2 with merons. This definition of merons is consistent with the semiclassical 
picture. In fact, an instanton contains two Wolff clusters, both with topological 
charge 1/2. At 6 = the clusters are independent and hence the merons form an 
ideal gas. The resulting disorder is responsible for the mass-gap of the theory. At 
6 = 7r the Boltzmann factor exp(iOQ) is 1 for even Q and —1 for odd Q. When a 
configuration contains a half-odd-integer charged cluster, flipping it changes Q by 
an odd integer. The resulting configurations then have opposite Boltzmann weights 
and their contributions cancel in the partition function. Therefore, at 6 = tt only 
those configurations for which all clusters carry integer charges contribute to Z{6). 
This means that the merons are now bound to pairs and can no longer disorder the 
system. Hence it is natural to expect that the mass-gap vanishes at 6 = tt. This 
agrees with Haldane's conjecture and it provides an exact formulation of Affleck's 
dynamical picture. 

We have performed numerical simulations with the meron-cluster algorithm at 
g = oo for volumes V = 3L x L = 18 x 6, 24 x 8, 30 x 10 and 36 x 12. For each 
lattice we have performed 10 7 sweeps for measurements. The correlation length at 
6 = is about 2.8 lattice spacings. Fig.l shows the topological charge distribution 
p(Q) on the 36 X 12 lattice. Due to our reweighting technique we can easily generate 
a distribution that covers 25 orders of magnitude. Up to logarithmic corrections 
discussed below we find Xt( 7r ) ^ L an d Xm( 7r ) ^ L which indicates a second order 
phase transition. For a first order transition the susceptibilities would grow pro- 
portional to L 2 . In fact, at very strong coupling one expects a first order phase 
transition [11, 9]. Due to the constraint in our action of eq.(l) the strongest bare 
coupling we can consider (g = oo) turns out to be too weak to reach this regime. 
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Figure 1: The topological charge distribution p(Q) on the 36 X 12 lattice. 

We have not run our algorithm at g < although this may be feasible. It would 
be interesting if one could then reach the hrst order domain. Affleck et al. have 
conjectured that the critical theory at the second order phase transition is up to 
additional marginal operators a conformal held theory in the universality class of 
the k = 1 Wess-Zumino-Novikov-Witten model [3, 9]. Close to the phase transition 
the mass-gap of the infinite volume theory should behave as 

m(9) = |#-7r| 2/3 |ln(|#-7r|)r 1/2 . (9) 

We consider Fisher's finite size scaling variable z = m(6)T which is a renormalization 
group invariant measure of the physical volume. Using the critical exponents of the 
Wess-Zumino-Novikov-Witten model with logarithmic corrections to scaling due to 
the additional marginal operators [3] one expects 

Xt (e,L) = T (lnZ)" 1 / 2 g t (z), X m{0,L) = T {InT) 1 ' 2 g m (z), (fO) 

close to the phase transition. Here gt(z) and g m (z) are universal functions. In figs. 2, 
3, Xt{@i £)£ _1 (ln L) 1 ! 2 and % m (^, Z)iv~ 1 (ln Z) -1 / 2 are shown as functions of z. In 
both cases the data are close to a universal curve. The logarithmic corrections to 
scaling are important for this. Note that no fitting or free parameters are involved. 
This confirms the described scenario. Cluster diagnostics on the 36 X 12 lattice at 
g = oo shows that most clusters are neutral, 4 percent have charge ±1/2, 1 permille 
have charge ±1 and very few have larger charges. The average sizes of the 0, ±1/2 
and ±1 charged clusters are 1.5, 10 and 32 lattice sites, respectively. 

In conclusion we have identified a second order phase transition at 6 = tt with the 
critical exponents of the k = 1 Wess-Zumino-Novikov-Witten model. This confirms 
Haldane's conjecture in the framework of the 0(3)-model. Identifying Wolff clusters 
of charge ±1/2 with merons we have provided a precise formulation of the meron 
picture. We have confirmed that the merons form integer charged pairs at the phase 
transition. The 0(3)-model is identical with the CP(f)-model. All CP( A)-models 
have instantons and hence a ^-vacuum. In the large N limit the CP( A)-model has 
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Figure 2: Data for the universal function gt(z). 
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Figure 3: Data for the universal function g m (z). 

a first order phase transition at 6 = tt [11] where CP is spontaneously broken. One 
expects that this persists down to N = 2 [11, 9, 12] although a numerical study in 
the CP(3)-model (motivated by the strong CP-problem) seems to contradict this 
[13]. Unfortunately, cluster algorithms do not work well for CP( iV)-models with 
N > 2 [14]. Lattice CP( iV)-models with N > 2 are of special interest because 
their topological effects have a well defined continuum limit g — > 0. This is not the 
case for the 0(3)-model [15]. Fortunately, in the condensed matter context of our 
investigation there is no need for scaling in the limit g — > 0. This limit corresponds 
to the classical limit S — > oo of the quantum spin system. The phase transition in 
6 is expected to persist for arbitrarily large spins, but there is no reason to expect 
scaling in S. 

Unfortunately the #-vacua in gauge theories cannot be simulated reliably at 
large 6 because sufficiently powerful algorithms are not available. However, the 
combination of cluster and reweighting techniques is applicable to the Ising model. 
This is presently under investigation. 




V = 18*6 

V = 24*8 
V = 30*10 
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We are indebted to I. Affleck and M. Liischer for very instructive communica- 
tions. 
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